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We report a first-principle theoretical study of the adiabatic decoherence undergone by a nuclear 
spin system in a solid, coupled to the phonon field through the dipolar interaction. The calculations 
are performed for a chain of weakly interacting 1/2-spin pairs, considered as an open quantum 
system in contact with a bosonic heat bath. By incorporating to the whole system Hamiltonian the 
fluctuations of the local dipolar energy produced by low frequency phonons, and assuming that this 
low energy fluctuations are adiabatic, we find that the spin dynamics can be described in closed 
form through a spin-boson model. The obtained results show that the coupling with the phonons 
destroy the spin coherence, and the efficiency of the process significantly depends on the complexity 
of the involved spin states. By using realistic values for the various parameters of the model, we 
conclude that this mechanism can be particularly efficient to degrade multi-spin coherences, when 
the number of ‘active’ spins involved in a given coherence is high. In this way, we show that the 
spin coherence in the adiabatic regime can be noticeably affected by this mechanism. 


I. INTRODUCTION 

Quantum dynamics of dipole interacting spin ensem¬ 
bles in solids arouses great interest in various fields of 
modern physics, both fundamental and applied. Partic¬ 
ularly, nuclear spins are suitable model systems to face 
frontier problems on the physics of many-body systems, 
like the very foundations of Statistical Mechanics [1]. 
Decoherence and irreversibility are essential pieces for 
the understanding of the complex dynamics which pre¬ 
cedes equilibrium, and they are ultimately linked with 
basic open questions such as the emergence of thermody¬ 
namic equilibrium from the underlying microscopic uni¬ 
tary quantum dynamics HE]. Within this context, is the 
challenging problem of explaining the mechanism that 
enables spins in solids to attain a quasi-equilibrium state 
over an early time scale -in its transit to equilibrium- long 
before the process governing thermalization may have 
acted. 

In the field of applications, considerable effort is ded¬ 
icated to manipulating and using quantum spin systems 
for new applications, like quantum computations and 
quantum information processing. A common character¬ 
istic of these developments is that they all need the oc¬ 
currence of quantum coherence between different states 
of a many-spin system |4|. Thus, a better understand¬ 
ing of the sources of the environment-induced destruc¬ 
tion of coherent superposition states, that is, of decoher¬ 
ence, becomes fundamental both in the search of scalable 
quantum devices and in the characterization of complex 
(many-body) quantum systems. 

A very sensitive technique to probe the nuclear spin 
dynamics along a wide range of well differentiated 
timescales is Nuclear Magnetic Resonance (NMR). To 
our knowledge, no treatment of the irreversible decoher¬ 
ence of solids in NMR has been reported that considers 
the interacting spins and the environment as a compos¬ 


ite system evolving under purely quantum rules. Both 
basic and applied research fields, would benefit from the 
development of Hamiltonian models which capture the 
essential physics and are amenable to rigorous analysis of 
the many-body nature of the quantum interacting parti¬ 
cles coupled to a quantum environment. To advance in 
this direction it is necessary to adopt a realistic model 
for the solid, which allows both to derive a detailed the¬ 
oretical description from first principles and to elaborate 
a quantitative analysis of the results. A theoretical ap¬ 
proach based on the usual spin-boson model 13 Ej that 
takes account of the spin interactions could meet those 
expectations. Particularly, this strategy could serve to 
shed light on the connection between decoherence and 
system-bath entanglement of a system of quantum inter¬ 
acting particles PEI- 

In this work we study the adiabatic decoherence of a 
system of dipole coupled spins interacting with a phonon 
bath. The term “adiabatic” means that spins and en¬ 
vironment do not exchange energy in the average. We 
follow the basic formal guidelines of the well-known spin- 
boson model [5|, however, instead of considering a set of 
uncoupled spins interacting individually with the boson 
field P [in] , in our model the observed system is a net¬ 
work of weakly interacting spin pairs, the boson bath 
corresponds to the phonons of the lattice at temperature 
T, and the system-environment interaction is generated 
by the variation of the dipole-dipole energy due to cor¬ 
related shifts of the spin positions from the equilibrium 
ones produced by the lattice phonons. For the sake of 
simplicity of the calculations we keep the main part of 
this interaction energy coming from the variations of the 
local (intrapair) energy. In other words, we consider an 
ensemble of spin pairs which are magnetically coupled 
with all the other pairs and are also correlated with each 
other through their interaction with the collective modes 
of the boson bath. In order to facilitate the calculations 
while keeping the main features of a real system, we con- 
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sidered a unidimensional chain. However, the extension 
of the analysis to other spin distributions is straighfor- 
ward. The strategy is based on calculating the exact 
quantum dynamics of an initial state of the observed sys¬ 
tem to derive the adiabatic decoherence function of this 
model. In order to compare with the experiment we es¬ 
timate the characteristic decoherence time of a typical 
element of the reduced density matrix. 

Our calculation aims to inquire if this pair-boson cou¬ 
pling can act as an effective source of decoherence, able 
to adiabatically bring an out-of-equilibrium system of 
weakly coupled pairs to a state of “internal” quasi¬ 
equilibrium. The occurrence of diagonal states represent¬ 
ing quasi-equilibrium states in solid state NMR mHH) 
is often postulated on the basis of the spin temperature 
hypothesis [IMS]- The spin dynamics associated with 
these states is described in the framework of ‘spin ther¬ 
modynamics’, by applying the tools of statistical mechan¬ 
ics to the spins in the quasi-equilibrium state. 

Formally, the statistical properties of the spin system 
is ‘kinematically’ described by extending to the quasi¬ 
equilibrium the mechanical statistical techniques. Thus, 
a Boltzmann operator involving the various (quasi-) con¬ 
stants of motion or quasi-invariants is postulated, as 
an extension of the equilibrium grandcanonical ensem¬ 
ble technique. However, the theoretical explanation of 
the transient processes occurring in the way of the sys¬ 
tem towards quasi-equilibrium still represents a challenge 
for quantum theories of many-body interacting systems. 
Explaining the decay of the off-diagonal density matrix 
elements could provide additional insight on the idea be¬ 
hind the spin temperature assumption. Similarly, this 
would allow to shed light on the influence of the environ¬ 
mental degrees of freedom on the generation of multi-spin 
correlation in the spin dynamics of the solid state . 

Section [ll| contains a brief review of the main general 
procedures of the theory of open quantum systems that 
will be used in the following sections. In Section [m] the 
spin system is defined, together with the Hamiltonians of 


the system, bath and spin-bath interaction. Section IV 


contains a derivation of the decoherence function and an 
estimation of a characteristic decoherence time scale. It 
is worth to stress that the simple model system (a chain 
of dipole coupled spins) used here is not aimed to charac¬ 
terize an actual crystalline sample; it is instead intended 
to serve as a means to prove that phonons can serve 
as mediators of an irreversible coherence loss. Finally, 
the Appendix A is dedicated to validate the adiabatic 
assumptions used in the calculation of the decoherence 
function. 


II. OPEN QUANTUM SYSTEMS 

As usual when describing an open quantum system in 
contact with an external bath, we write the total Hamil¬ 
tonian as 

( 1 ) 


where the system of interest is described (in isolation) by 
Hs, the bath hy Hb and Hsl represents the system-bath 
interaction. 

Since we are interested in describing processes which 
take place without energy exchange between the spins 
and the bath, the Hamiltonians must satisfy the adiabatic 
condition 


[Hs,Hsi]=0 (2) 

which means that the mean value of the spin energy, {Hs) 
is a conserved quantity. This condition is essential in this 
approach since it rules out relaxation and thermalization 
effects. 

To model the interaction Hamiltonian we follow ref. 
P5] and write 


USL 

A 

as the sum of products of generic Hermitian operators, 
A"^', of the observed system (spins) and which acts 
on the Hilbert space of the bath. The only restriction 
imposed to Hsl is condition (H): [Hs,Hsl] = 0. 

The properties of the observed system are expressed 
by a reduced density operator a, which is obtained by 
tracing the density operator p of the whole system over 
the bath variables [24], that is, 

cr = Tib{p}- 

The strategy for studying decoherence involves calculat¬ 
ing the time dependence of a. The unitary dynamics 
of p is driven by the time evolution operator U{t) = 
exp{—ftH}, which, owing to condition ([^ can be fac¬ 
torized as 

U{t) = (3) 

where we introduced a “rigid lattice” time evolution op¬ 
erator Vo{t) = which acts on the spin variables 

only and an operator V {t) which acts on both the spin 
and lattice variables. Let us further define an operator 
Vm{t) acting on the lattice variables only as 

V{t)\m) = \ m)Vmit), 

where | m) is an eigenvector of Hsl- 

The time dependence of an element of the reduced den¬ 
sity matrix is 

cTmnit) = {m\TiB |c/(f)p(0) I n). (4) 

We now make the usual assumption that the initial con¬ 
dition can be factorized as p(0) = cr(0) 0 Pb, then, 

(Tmn{t) = {m\Voit)a{0)V^{t)\n) Ttb {Vm{t)pB V^{t)} 

( 5 ) 

This very general expression may now be applied to 
particular cases by selecting the operators 'Hs,Hsl and 


H = Hs + Hb + Hsl, 
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FIG. 1: This figure represents the model used in the calcula¬ 
tions. We assume a monoatomic chain with lattice parameter 
a, with two atoms per unit cell, separated a distance d. The 
chain is embedded in a magnetic field which form an angle 9 
with the chain. Each spin interacts with all other spins of the 
sample through dipolar interaction. 


'Hb- The factor which involves ct( 0) in Eq.(^ describes 
the complex (unitary) dynamics of a closed interacting 
system, while the non-unitary behaviour of this time evo¬ 
lution is given by the trace over the bath variables. Our 
aim is to calculate the second factor in the case of a linear 
chain of dipole coupled spin 1/2 pairs interacting with a 


bath of phonons, as described in Section III 


III. MODEL SYSTEM 


We consider a system of dipole coupled spin 1/2 pairs 
in a strong, static magnetic field and represent the model 
solid as a monoatomic linear chain with a basis. The 
primitive cell of this unidimensional Bravais lattice with 
parameter a, contains two identical atoms of mass m at 
a distance d < a/2. There are two elastic constants K 
and G between neighbouring atoms (with K > G). See 
Fig. 0. 

The phonon bath describes the small amplitude dis¬ 
placement of nuclei around their equilibrium positions in 
the crystal which in turn perturbs their magnetic dipole- 
dipole interaction. The bath energy of this chain model 
(in units of h) can be written as the sum of uncoupled 
oscillators of frequency ujk (neglecting the zero point en¬ 
ergy) as 

Ub ='^tOkblbk, ( 6 ) 

k 


where the sum runs over the first Brillouin zone and b ^, b/ 
are the annihilation and creation operators, respectively, 


which satisfy 


bh, b 


= 1. The dispersion relation of a 


ID chain [25] is 


± + G^ + 2KGcosika), (7) 


which admits an optical and an acoustic branch. 

In this model each spin interacts with both an external 
magnetic field (Zeeman energy) and with the other spins 


of the sample through dipolar interaction. The Zeeman 
term is 


= = + If ), (8) 

A A 

where A labels the pair in the chain, and if is the 2 - 
component of the angular moment operator of the m— th 
spin in the pair A. 

As regards the dipolar term, He, we first adopt the 
high magnetic field approximation, since we assume that 
1177211;^ WHdW and therefore keep only the secular part of 
the dipole Hamiltonian, which commutes with Hz (this 
assumption is very usual in NMR). Owing to the geom¬ 
etry (chain of pairs) selected in this work, a hierarchy of 
dipolar couplings arises which naturally allows to split 
the secular dipolar Hamiltonian into two terms, one in¬ 
volving the interaction between spins belonging to the 
same cell, which we call intrapair, and the other repre¬ 
senting the mterpair interactions, 710= + 

Then we write 

^^ntra = ^ ^ ^ ^ {sif if - ( 9 ) 

A=1 A=1 


with defined as 

S2D = -5^[3™y«)-l], (10) 

where 9 is the angle between the external field and the 
chain direction, yp the proton gyromagnetic ratio and 
rA is the distance between the interacting nuclei of the 
A pair. In the following calculations we will not need 
the explicit form of the interpair term, however, 

we assume that it can be truncated so that it commutes 
with the intrapair term |26| . 

We have so far defined the spin interactions and the 
bath energy. The connection between these two subsys¬ 
tems is provided by the change in the dipolar energy gen¬ 
erated by the lattice vibration. It is worth to notice that 
it will only affect the dipolar and not the Zeeman inter¬ 
action. We now assume that this perturbation is small, 
so that the dipolar energy can be expanded around the 
equilibrium value of the internuclear distances. In the 
case of the intra-pair interaction, we have 

(rA-d). (II) 

rA=d 

The first term of Eq. 0 is the rigid lattice value of the 
dipolar intrapair Hamiltonian of pair A. The second term 
is the correction due to lattice vibration. Since the aim 
of this work is to examine if the pair-boson interaction 
may in fact be a source of decoherence, at this stage, 
and for the sake of simplicity, we consider only intra-pair 
variations. This assumption is also based on the fact that 
■^mter jg gjj^aller than because of the dependence 


nj A 

Hd 


n 


dAf 

D\rA=d~^ 
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r ^ (see Appendix A for more details). That is, 


'Hd 


[njintra , njinterl , \ 

[HOg + J+2^ 

A 


dU^ 

dr A 


{rA-d), (12) 

rA=d 


where Hdo corresponds to the rigid lattice. 

We write {ta — d) in terms of the displacements from 
equilibrium, ^a,u, of each spin u = 1,2 at molecule A, 
since they can be easily related to the phonon creation 
and annihilation operators bk^s and b\. ^ of mode k and 
branch s of an A^ primitive cells chain, as |25j 



as) 

where (A ■ a + u ■ d) is the equilibrium position of spin u 
in pair A. Then 


rA-d = HA,2 - MA.i = {9k*bk + gk^l) , (14) 

k 


where we omitted the branch index s to simplify notation, 
and defined the coupling constants 


^—ikAa^ 

9 k — 6 9k^ 


(15) 


with 


= (1 - 




h 

2muJkN 


In this way the term in square brackets in Eq.(12) is 
evaluated at the positions of the rigid lattice and the sec¬ 
ond term, having both spin and lattice variables, emerges 
as the system-bath interaction. We then define 


ns ^Hz + EAn 


D\d 


njinter I 


nsL =Ea^^^ Ek {g^bk + 


where 


A^ = 


dH^ 


dr 


rA=d 


Hd 

'~d 


A1 tA2 _ jAl jA2 




(16) 


(17) 


is an operator acting on the spin variables only, and Hd = 
Hd l»'A=d' 

Finally, we are going to assume that = 0, 

wich is equivalent to consider low-frequency modes in the 
interaction Hamiltonian (see Appendix A for detailed dis¬ 
cussion). This allows to factorize the evolution operator 
as in Eq. (|^. 

It is worth to mention that the Hamiltonian Us in¬ 
volves the static magnetic dipole interactions between all 
the spins, both intra-pair and inter-pair. In this way, the 
time evolution under Us only, still reflects the complex¬ 
ity of a multiply connected dipolar network. It is only in 
nsL where we neglected the variation of the inter-pair in¬ 
teraction with the displacements. As a consequence, the 
variations associated to which are not considered 

in this work, would add an extra source of decoherence. 


IV. DYNAMICS 


In this section we calculate the exact coherence dy¬ 
namics given by Eq.([^, using the Hamiltonians of the 
particular system defined in Eqs.Q and We also 

assume that the initial state of the composite system is 
separable and that the bath is in a thermal state, de¬ 
scribed by independent density matrices 0^ 


p(0) = (t( 0) (g)]^0fe 

k 


with 


Ok = 

Zfc = (1 - 


(18) 


(19) 


In accordance with our definition of nsLt there exists 
an eigenbasis {|m.)} whose elements can be written as 
direct products of states of the different pairs. 


|m) = \m)i<Si\m) 2 <^ = \ mi,m 2 , ■■■,171^), 

( 20 ) 

where each Iuia} is any of the four states of pair A, that 
is, eigenstates of Hd and 

|mA)e{|l,l),|l,0),|l,-l),|0,0)}. (21) 

The corresponding eigenvalues of A"^ are 

A^Itoa) = Arn^lm^). (22) 


The action of Hsl over a state \m) is 


'Hsl m) 

= Ea^^® Ek (g^bk + gtbl 

j |mi, ..^mN) 


~ EjA Am^ 

Ek + fi¬ 

^'‘^‘^gkbi) \m) 


= Ek 

fe 9k bk -f ^m,k gk ^fe 

) \m), 



(23) 

where we 

defined 




^m,k 

^J2e-^kAa^^A. 

(24) 


to emphasize that the sum over A in Eq.(23) links the 
eigenvalue of each pair with its position in the chain. 

We now introduce the notation 


nsL (m) = {Kn,k9kbk + > ( 25 ) 

k 


which allows writing 
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(remember that = \m)'HB)- From Eq.Q and 

Eq.(j2^, the exponent at Eq.(26) can be written as a 
sum over different modes k as 

T-Lb +'HsL{,m) = ^2lk,m, (27) 

k 

where 

')k,ra = ^k^ppk “t” ^m,k9k^k “t“ ^ni,k9k^\ 

Since terms with different k commute, the time depen¬ 
dence of the density operator of Eq. ([^ becomes 


CTmnit) = {m\Vo{t)(j{0)Vfj{t)\n) 

xrifelTrfc 


(29) 


where the trace over the phonon modes k comes from 
the trace over the states of the heat bath. Equation 
(29) has the same structure of Eq.(2.17) from Ref. [TU] . 
but here the coefficients ^k,m have a definite meaning in 
terms of the lattice parameters. The key problem now, 
is how to calculate such trace. Some authors solve spin- 
boson problems by calculating an explicit expression for 
the time evolution operator 13 [27], an alternative strat¬ 
egy uses the coherent states for the harmonic oscillator 
to calculate the traces US HE). We adopted the latter 
strategy. 

Coherent states | z) are eigenstates of the annihilation 
operator b^. The trace of an operator in the basis {| z)} 
can be calculated with the following integral [55|, 


with 


Using that 


Tr(0) = J d^z{z\0\z), 


d^z = —d(Re z) d(Im z). 

TT 


(30) 


(31) 


1 = jd^z\z){zl 
we get for arbitrary k 

Trfc 0 = 

= ^/d2zid2z2(i2z3(zi|e-*T'''.-*|z2) (32) 

x(z2|e-^^‘^'''’fc^''|z3)(z3|e*T''=.-*|zi). 

It is well known that [55] 

(Z2|e-^^‘^''''^'"=|Z3) = (33) 

In order to calculate the two other matrix elements of 


the complex exponential operators in Eq.(]32[), it is con- 

irs 

r] = b+'-^ (34) 


venient to introduce the Bogoliubov shifted operators 

^m9 


which allows writing 


t \>^m?\9? 

7m =ujri^r] - 

UJ 


to obtain 


(zi|e-*T'’-‘|z2) = e*‘ 




{Zl\z2) 




(35) 


(36) 


In this way, calculation of the trace in Eq.(32| implies 


calculating six gaussian integrals. The result is analogous 
to the one reported by Ref. m, that is, 

kmn(f)| = \{m\Vo{t)a{0)V^ (37) 

where Tmnit) is the decoherence function 


(38) 


Tmnit) = ‘2j2ks\^rn,k,s- K,k,s\^ 

’ ^k,s 

X (1 — COSUJk,st) coth 
(we restored the branch index s). 

A. Analysis of the decoherence function 


In this section we work in the expression of the deco¬ 
herence function in order to analyze its time dependence 
and estimate a decoherence time in terms of realistic val¬ 
ues of the lattice parameters. 

In the continuum limit we integrate over k > 0 due to 
the parity of the integrand, so, Eq.(38) transforms into 


Fmn(0 ~ d / ^ '^n,fc 

„ JO 


|2 \9{^^k,s)\'^ 


kOl 


(39) 


where 


(1 — cos ijJk,st) coth ^^^'^^ G{k)dk 


(1N 

G{k)dk = —dk 

TT 


is the density of states of a linear chain and 


\9M\^ = 


mNuik,: 


(1 — cos kd). 


(40) 


(41) 


Since the integral in Eq. (39) cannot be calculated analit- 


icaly, we now analize the behaviour of the integrand as a 
function of k. In first place, we consider only the acoustic 
modes since they are the most populated at room tem¬ 
perature, then we use the dispersion relation 


KG 


UJk 


2m{K + G) 


ak = ck 


(42) 
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and write the integrand as 

coth (^ (43) 

Y 2 y mn ck 

= LlJik)hik,t). 

We defined the functions Lmru h{k,t) and f{k) to con¬ 
sider their behaviour separately 


• The dependence of 


k-'mn — |A77T,^fc A^^^/el — 


E' 


— ikAa 


{^niA A„^) 


on the spin states | m) and | n) is complex. A di¬ 
rect calculation of L^n would involve a sum over 
the 4-^ states, then, we must seek for convenient 
approximations to find an estimation of its possi¬ 
ble values. 

The four eigenvalues A™^ of operator (from 
Eq.(17)) of each A pair are 


A^|l,l) =0.5^11,1) 
A^|1,0) =-lf^|l,0) 
A^|l,-1) =0.5^|1,-1) 
A^|0,0) =0 


(44) 


From Eq.(39) it is clear that Emm = 0 and conse¬ 
quently Tmm = 0, which is a consequence of ([^. It 
is also possible that Lm,n = 0 for m ^ n for some 
particular states, like | to) = |1,1) 0 |1,1)...|1,1) 
and I n) = |1, —1) 0 |1, —1)...| 1, —1). In a more 
general case may involve up to N nonzero 

terms. Let us assume for example that only M of 
the terms contribute to the sum over A and also 
that the coefficients (Am^ — A„^) ^ In such a 
case 


7-2 


„—ikAa 


E ivi 

A^ 


2 l — cos[Mka] 
\ d J 1 —cos(fea) 


(45) 


h(fc,t) = = ^sinc^ (f*) 

takes its maximum at k = 0 and its first 

2tt 

minimum (as a function of k) at fco = —• 

ct 

Concerning f{k) = coth ^(1 - cos(fcd))^, 

there exists an interval (O^kmax) where 

-1 


coth^ - (^) 
I — cos{kd) ~ 


(46) 


is valid, and the function /(fc) can be approximated 
by its small-fc expression, f{k). 

We now define 

i{k,t) = Lmn f{k) h{k,t), 

which should approximate I{k,t) provided kmax ^ 
(because h{k,t) becomes very small if fc 3> ko). In order 
to analyze the time range where this approximation is 
valid, we compare I{k,t) and I{k,t) at different times. 
To do so, it is necessary to set the values of the different 
constants involved. We choose d = O.Inm and a = 3d, 
which are the usual values for interatomic distances in 
solids (for example the distance between Hidrogen nuclei 
in a hydration water molecule); to = 1.66- I0“^^kg, is the 
mass of an Hidrogen atom, = 300kHz and T = 300K 
and for the speed of sound in a solid we use c = 3000m/s. 

We can now proceed to calculate the exact and the 
approximate integrands, assuming, for example, M = 
10^. Figurej^shows the plots of I(k, t) ^lid) and /(fc, t) 
(dashed), calculated at t = Ips in Figj^a) and t = Ins 
in Figj^b). It is evident that the two functions coincide 
in the lower plot, showing that the condition kmax ^ 
is valid for times t > Ins. 

Therefore, under these conditions it is possible to ex¬ 


tend the upper limit of the integral in Eq. (39) to get an 
approximate decoherence function 


r = 


2rtjjkBTaM'^ f°° 1 — cos(cfct) 


nine* 




dk 




t 


(47) 


mc^ 


~ 2.5 X t. 


This expression exposes the fact that Tmn depends 
on the number M of spins involved in the transition 
which links | to ) and | n ) states. 

In the limit of small k, 


limL^„ = L^„=l^) 


This numerical result implies that coherences in this par¬ 
ticular spin system, decay in a characteristic time within 
the millisecond range if M = 10^, which is an experi¬ 
mentally accessible time scale. Besides, a most relevant 
feature of our result is that the decoherence rate is pa¬ 
rameterized by the eigenvalues of the spin part of the in¬ 
teraction Hamiltonian. Eq.(47) contains also an explicit 
dependence on temperature, however, there can also be 
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integrand by its low k value, I{k,t) which results in 
1 — cos(cfct) 




fc2 


-dk 


— 1 + cos(ect) + Si(ect) ect 


(50) 


where Si stands for the sine integral. Since the second 
line of Eq.(50) tends to infinity when t —>■ oo, then 


^ , , —1 + cos(ect) + Si(ect) ect 

Tmn{t) > -^^^-^ 00 


if t —)■ oo 
(51) 



FIG. 2: (color online) Plot of function I{k,t) from Eq.(43l 
(solid line) and the approximate function I{k,t) (dashed) for 
two different times; (a) t = 1 x 10~^^s and (b) t = 1 x 10“®s. 
Both lines coincide in (b). 


an implicit dependence of other parameters for example 
due to the finite phonon lifetime and the speed of sound 
in the particular system. 

Although we made a realistic approximation, it is im¬ 
portant to show that the decoherence function diverges, 
independently of such approximations and of the chosen 
values of the parameters. 

Lemma.. lim^_^QQ 7^ 0- 

Proof Let us first split the contributions from acoustic 
and optical phonons. 


r = r 

— -L n 


r 


op 


(48) 


Since both terms are positive (because their integrands 
are positive), it suffices to prove that any of them diverge, 
then we consider Foe. In order to find a lower bound 


for the integral of Eq.(39l, we see that for an arbitrary 
wavelength e. 


t/q 


"acO: / {...)dk- 

Jo 


{...)dk> / (...)*. (49) 

Jo 


By setting e small enough, it is possible to replace the 


This interesting result shows that the matrix elements 
of Eq.(29) attenuate because the lower bound for r{t) 
attains large values as time grows, so proving that de¬ 
coherence due to the interaction of a chain of dipole- 
coupled-pairs of spins with a phonon bath is possible. 


V. DISCUSSION 

The main result of this work is the finding that a chain 
of weakly coupled spin pairs, adiabatically interacting 
with a phonon bath, treated as an open quantum sys¬ 
tem, can indeed undergo irreversible decoherence over 
an intermediate time scale, earlier than that dictated by 
thermal processes. The possibility of completing (attain¬ 
ing) a closed derivation of the decoherence process in a 
realistic model gives us the opportunity of highlighting 
the main features of this subtle topic that can be ex¬ 
pected to find in a real system. Therefore, this interest¬ 
ing academic example can provide insight into the prob¬ 
lem of irreversible spin dynamics in solids. The proposed 
model can be generalized to be applied in a real system 
of weakly coupled pairs, like hydration water molecules 
in hydrated salts, and could also be extended to ordi¬ 
nary solids. The explanation of irreversible decoherence 
in solids has remained as an open question for a long 
time in the field of solid state NMR. Particularly, the 
mechanism by which nuclear spins of a solid are able to 
achieve a state of quasi-equilibrium continues to be elu¬ 
sive nowadays. The conclusions obtained in this work 
can shed some light on this question. 

The formal aspect of our work is inspired by the spin- 
boson model strategy for treating the time evolution op¬ 
erator in the case where the interaction and bath Hamil¬ 
tonians do not commute. In this framework it is possible 
to derive the time dependence of the reduced density ma¬ 
trix in ‘exact way’, without resorting to coarse-graining 
procedures. The traditional spin-boson model consid¬ 
ers individual spins interacting with a boson field, and 
the spin-environment interaction is represented through 
a quantum Hamiltonian linear in the boson variables, in¬ 
volving phenomenological coupling strengths, which nat¬ 
urally satisfies the adiabaticity condition given the form 


























of the system Hamiltonian. The effect of the unkown 
coefficients on the result is introduced a posteriori by 
modeling the spin-bath interaction. In contrast, in our 
proposal the system is thought as a network of dipole- 
dipole interacting spins, whose coupling with the envi¬ 
ronment originates in the small variations of the local 
spin interaction produced by lattice phonons. Thus, the 
system-environment interaction is represented by a first- 
principle dipole-phonon Hamiltonian, in which the physi¬ 
cal magnitudes are expressed in a fundamental way. This 
characteristic of the model may facilitate a close com¬ 
parison with experiments. Even when the variations 
of the inter-pair dipolar interactions were not included 
in the system-environment Hamiltonian, the mechani¬ 
cal correlation of the different pairs established by the 
phonon displacement field is reflected in the final re¬ 
sult as a collective feature of the decoherence process. 
A similar discussion was given in the study of decoher¬ 
ence in nematic liquid crystals , where even when the 
molecules are magnetically isolated in average and the 
main contribution is calculated using the intramolecular 
spin-environment Hamiltonian, the calculated decoher¬ 
ence function reflects the influence on a given molecule 
exerted by the other molecules through their mechanical 
interaction as a quantum effect |301131j . 


It can be seen from Eqs. (§ and ( [ 3 ^ that the deco¬ 
herence mechanism and the multi-spin dynamics of the 
observed system are complementary aspects of the irre¬ 
versible dynamics of a system observable. Because of 
the multi-spin character of 'Hs, states |m) and \n) are 
not eigenstates of the evolution operator Vo(i)j then the 


closed-system dynamics represented by Vb(i)o'(0)lo^(t), 
connects eigenstates of Hj where the number of “active” 
spins (those which change their state in an m —>■ n tran¬ 
sition) can be large. This correlation grows with time 
[20l [32] at a rate which depends on the particular dipolar 
network. 


By other hand, the decoherence function Tmn also de¬ 
pends on the number of active spins and is responsible of 
the irreversible decay of the matrix element amn- This 
microscopic irreversibility could convey an irreversible 
behaviour to an observable {0(t)) = Tr [a{t)0]. In NMR 
experiments the observed signal decay during the evolu¬ 
tion under the spin interactions is usually interpreted as 
a consequence of the transformation of single-quantum 
coherence terms to unobservable multiple-spin, single¬ 
quantum coherence terms, during the unitary evolution 
of the initial spin density matrix in a closed system |33j . 
In such case, the system evolution would be completely 
reversible and it would be possible, in principle, to per¬ 
fectly refocuse the spin dynamics. Due to the same 
reason, a quasi-equilibrium state would not be possible. 
However, the system is not completely isolated, and de¬ 
pending on the efficiency of both processes, namely the 
correlation growth and the decoherence effects, the mi¬ 
croscopic information leakage can produce an irreversible 
macroscopic signal decay. The estimation of the decoher¬ 
ence time scale was made by using realistic values of the 


various parameters involved. The order of magnitude 
obtained in Eq. (471, showed to be consistent with deco¬ 
herence rates measured in NMR refocusing experiments, 
like the ‘magic echo’ and MREV pulse sequences [2011311 ■ 
Also, the obtained values are compatible with the time 
scale associated with the generation of quasi-equilibrium 


This suggests that adiabatic decoherence due to the 
pair-phonon interaction should also be considerded as a 
source of authentic irreversible attenuation of refocused 
echoes in NMR experiments. This result contrasts with 
the null expected influence of dipole-phonon interaction 
in the processes of spin relaxation and thermalization 
[I71|31|3S|. This fact should not be surprising since 
decoherence and thermalization are processes of very dif¬ 
ferent nature. Thermalization is generally described by 
Markovian master equations for the reduced density ma¬ 
trix that evolves in a coarse-grained time scale, with a 
rate controlled by thermal fluctuations of the spin-lattice 
interaction. On the contrary, adiabatic decoherence is a 
full-quantum process which involves, as a prime ingre¬ 
dient, the correlation between system and environment 
generated by the non-commutation of the interaction and 
bath Hamiltonians, ['Hslj'Hb] ^ 0. 


We saw that the low frequency phonons are which play 
the main role in the decoherence function. Besides, the 
spin-bath coupling associated with these modes guaran¬ 
tees the existence of an adiabatic time scale where the co¬ 
herence is lost without appreciable influence of thermal 
fluctuations. This strong relation between decoherence 
and low-energy excitations remains to be investigated in 
solid state NMR. In summary, we found that the dipole- 
phonon interaction provides an effective mechanism of 
adiabatic decoherence for a system of weakly interacting 
spin pairs in the solid state, even when the thermaliza¬ 
tion times can be very long. The model presented allows 
the estimation of decoherence rates in terms of parame¬ 
ters of the system. The system-environment coupling is 
introduced through the variations that the phonon field 
produces on the intra-pair interaction energy, considered 
as the main contribution. Though the interaction Hamil¬ 
tonian is a sum of individual terms, the strong correlation 
that exists between the pairs due to the phonon field is 
a decisive feature of the decoherence function (like in the 
case of liquid crystals). 


The adiabatic decoherence induced by the coupling 
with the environment can be considered the initial ir¬ 
reversible mechanism responsible for the degradation of 
the signal in the intermediate time scale, long before the 
thermal fluctuations lead the system towards the final 
equilibrium with the bath. This mechanism can explain 
the signal decay in refocusing experiments and also the 
occurrence of quasi-equilibrium states over the interme¬ 
diate time scale. 
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Appendix A: The adiabatic approach 


This appendix is dedicated to validate the adiabatic 
approach used in the calculation presented. In such ap¬ 
proach, the Hamiltonian of the dynamics can be replaced 
by an adiabatic one, which rules the dynamics of the sys¬ 
tem in the decoherence time-scale (i.e. the time-scale far 
earlier than the thermalization one). 

The section is divided in four parts. In |A ![ we present 
the definition of the complete Hamiltonian of the dynam¬ 
ics and the chosen writing for the parts that compound 


it. In A 2 is described an analytical calculus of the ap¬ 
proximations used for the dynamics. In |A 3[ we discuss 
who is the adiabatic Hamiltonian of the dynamics, which 
is obtained from the complete Hamiltonian. 


1. Description of the complete Hamiltonian of the 
dynamics 

The complete Hamiltonian of the dynamics is formed 
by the following Hamiltonians: Zeeman, complete dipo¬ 
lar, the interaction and the bath or environment one, that 
is 


n = nz + nD + ni + nB, ( ai ) 

which will be described hereupon. 

The Zeeman and the dipolar Hamiltonians applies only 
over the spin-space states, we can write them as 

(A2) 

nD=Yl (A3) 

n——2 

where 

(A4.) 

A 

E E <Sa..b^y 

A;u=l,2 B^A;v=l,2 


spins (or the system) and the bath (or the environment 
or lattice), respectively, and 1 is the identity operator. In 
Eqs. (|A3|) and (|A4|), the subscript Da and De represents 


the intra and the inter parts of the dipolar Hamiltonian, 
on the other hand, the superscript n indicate its compo¬ 
nent 0,±lor ± 2. The indexes A and B run over the 
spin-pair number, u and v run over the number of spin 
inside each pair (defined by the numbers 1 and 2 ). 

The interaction Hamiltonian is derived from the dipolar 
one and can be expressed as 


2 

ni= Y. + (A5) 

n——2 

where 

A 


njn — 

^7e= 2 


E E 

A-u-1,2 B^A\v=l,2 


n 


n{s) 

Ie{A,u;B,v) 




(&) 

Ie{A,u\B^v) ■ 


(A 6 b) 

In the same way that Eq. (A3), in Eqs. (A5) and (A 6 ) 
the subscript la and le represent the interaction part 
derived from the intra and the inter parts of the dipolar 
Hamiltonian, and the superscript n indicate the corre¬ 
sponding component. 

Finalizing the descriptions of the Hamiltonians, the bath 
Hamiltonian is writing as 


•Hs = (A7) 

Now, a detailed expression for the Hamiltonians previ¬ 
ously shown will be presented. The well-known Zeeman 
Hamiltonian is 


nf^u^YtiAY (A8) 

A 

where cug = jBq is the Larmor frequency, 7 is the proton 
gyromagnetic ratio and Bg is the strength of the static 
magnetic field, which is applied along the laboratory z 
axis, and 1 ) 2 ) ^ projection of 

the total proton spin angular momentum of the A-th 

(s') 

spin-pair, with as the corresponding angular 

momentum of the i-th spin inside each pair. 

To obtain the expressions of the different parts that 
conform the dipolar and interaction Hamiltonians, we 
will define several functions and operators, which are 
written under the international system of units (SI). 
First, we define the following functions 


Hereafter, the symbols of the form and indicate 
operators acting exclusively on the Hilbert space of the 


o,, /io7|>^ [1 - 3cos2(6»)] 

-- 


(A9a) 
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Ho^ph sin(20) 

Stt V 2 7-3 ’ 


(A9b) 


^±2{r) = - 


Mo7p^ 

Stt 


3 sin2(6») 
2 


(A9c) 


Then, the following tensors are defining 


rp(s) 

^ {2,0}(A,u-B,v) 


[3Iz(A,tt) Iz(i?,D) 


- f . f 1 


(AlOa) 


^Ie{A,u-B,v) 


k,l 


(A,u-B,v) .1(6)] 

yk,l '^k,l \ ’ 


(A15b) 


with 


iA^u\B^v) 

9k,I 


h 


2mujk iN 


^i2kd 


\y [^—'ik(Ba-\-vd) —ik(Aa-\-ud)' 

X c- — c- 


(A16) 


Finally, neglecting the zero point energy, for the bath 
Hamiltonian we have 


rp(s) 

^ {2,±l}(A,u-B,v) 


T 2 '^dz{A,u)^z{B,v) 


+ 1. 


iA,u)i-±(B, 


](^) 

u)J ) 


(AlOb) 


( 6 ) _ 


rt D — 


k,l 


bt(6) b(^) 
^k,l ^k,l- 


(A17) 


'^{2,±2}{A,u-,B,v) — 2 [^±(^.“)^±(S,i;)] ) (AlOc) 

also, we are using for the tensors that involve spins into 
a pair the following notation 


2. Leading evolution operator 

In this section an approximation of the evolution 
operator will be presented. Such approximation has the 
feature of being very close to the actual operator, under 
the condition that will be discussed. 


— rT.(s) 

{2,n}{A) — {2,n}(A,l;A,2)' 


(All) 


Therefore, using the previous definitions of functions 
and tensor, we have for the intra-pair dipolar Hamilto¬ 
nian 


First, we define a general Hamiltonian "H as the sum 
of two parts. One main part Ha, which has a norm or 
eigenvalues with similar order of magnitude that H, and 
other perturbative one T-Lp, which has a norm or eigen¬ 
values smaller than Ha and H. Therefore, we write 


(A12) 

On the other hand, for the inter-pair dipolar Hamilto¬ 
nian we have 


'^De{A,u-,B,v) — '^^^('''^•''^■•B,v)'T\2,n}{A,u;B,v)- 

In the same way, the parts of the interaction Hamil¬ 
tonian can be writing as follow. The intra-pair dipolar 
Hamiltonian have the following expression for the spin 
and the bath operators 


^la(A)- ^^Da{A)- ^ ^ {2,n}(A)’ 

(A14a) 


H 


(6) = 
Ia{A) — 


^ ^ ^—ikAa 
k,l 



(A14b) 


Otherwise, for the inter-pair dipolar Hamiltonian we 
have the following expression 


njn(s) 

^ I e{A^u\B ,v) 


3 

^ A,u-,B^v 


-T/"(s) 

De{A^u\B^v) ’ 


(A15a) 


H = Ha + Hp. (A18) 


Extracting the Hamiltonian Ha from the complete evo¬ 
lution operator, namely U{t), we can factorized as 


U{t) = e-^^'^ = e"** fo dr Hf, (r)] ^ 

where Hpir) = Hj) and e represents a de¬ 

velopment of an evolution operator with integrals which 
have been ordered in time, where the bigger the possible 
values of the time in the extreme of the integrals, more 
to the left the operator depending of that time will be. 
That is, such ordered evolution operator represents the 
following development 


= 1-i f dTiHfsin) 

Jo 

- f dTi f dT2Hli{Ti)Hp{T2) 

Jo Jo 

-\-i f dn f dT2 [ dT3H/3{Ti)Hi3{T2)Hp{T3) 

Jo Jo Jo 


+ i-^y 


dr-i ■ ■ ■ 


dTnHp{Ti) ■ --H^iTn) 


(A20) 



11 


Assuming that the Hamiltonian have a part that com¬ 
mutes with %a, namely and other one which does 
not do it, namely "H/j, then we can write 


we define an eigenbasis {| a)} of the operator T-La where 
'Ha\ct) = a\a). Accordingly, we calculate the value of 
integrals like the following 


'H.p = Tip + 'Hp. 


(A21) 


Using the definition (A21| in the integrals of Eq. (A20l, 
and defining 


= (A22) 

we will obtain for the first three integrals the following 
development 

dTiUpin) = npt (A23a) 



, a 


Uo 


dr fir) Hpir) 


\a') = 


(A25) 


[alnpla') / dTf{T)e’’' 


with 0Ja^a> = a—a'. It is worth to note that ( a {Hpl a) = 
0, because that Hp is the non-commutative part of Hp 
with respect to Ha ^ therefore always we have that a ^ o! 
in the integrals. 

If the exponential function varies faster than 

/(t), we can approximate 


fdn r dT 2 np{T,)np{T 2 ) = 

Jo Jo ^ 

+ {Hp t) J° it) + Xj it) (Hp t) - {Hp t) it) (A23b) 

dTl f dT2'HpiTl)'HpiT2), 

Jo 



Iit)= f ~ lim/(r) / 

Jo ^^0 Jo 

= tlim/(r)e*---‘/^Sinc(^^), 

(A26) 


[ dTl f dT2 [ dT3HpiTi)'HpiT2)'HpiT3) 

Jo Jo Jo 


iUpt)^ , iHpt) 


+ 


3! 

iUpt)^ 


2 X^it)-iHpt)^X^it) 
- iHpt)^ 


X^it) +X^it) 


2 > 2 
+ iUp t)Xlit) iHp t) - iHp t)Xlit) iUp t) 

Hiilpt) f dTl f dT2'HpiTi)'HpiT2) 

Jo Jo 


-iUpt) 


dTiTiHpiTi)Xj^iTi) 


+ - I dTlTll-LpiTl)i'Hpt)X^iTl) 


- dTinpiTi)i'Hpt)Xp (ri) 

Jo 

+ [ dTinpiTi)X^iTi) iHpt) 

Jo J 

+ [ dTl f dT2 [ dT3HpiTi)HpiT2)HpiT3), 

Jo Jo Jo 

(A23c) 

where is defined 


X^it)^-j^dTT-HpiT). (A24) 


where is defined 


Sine 


with the properties 


,a't\ _ siniu}a,a't/2) 

UJa,a'tl2 


Sinc(O) = 1, 


Sine 


O^a.a't 


~ 0 . 


(A27) 


(A28a) 


(A28b) 




For the cases of the integrals in Eq. (A23l, we have 
that lim 7 ._j.o /(t) = 0, in particular /(t) = (r/t)” 
(n = 0,1, 2,...) for (A24|. In the cases of our concern, 


we will have that Wq q,' ;g> 1, thus 




Using Eqs. (|A28 ) and (A29l, we obtain for the integrals 
in Eq. ( A23pthat 


r = e--<^'*/2tSinc(!^) , n = 0, 
1-0,717^0. 


The value of the matri x elements of the operators of 
the integrals in Eq. (A231 can be calculated and thus we 
have an order of magnitude of such integrals. In this way, 


Therefore, we conclude that it is possible to neglect 
for all t the values of the integrals that involves powers 
of r, with respect to the other one which does not do it. 
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Accordingly, we can very well approximate 


reference for the dynamics. 






3! 


dTiT-Lpiri) 


~t“ ( 1 — 'i'fi/3 t “t“ 


f dTi [ dT2'Hl3{Ti)'Hp{T2) 

10 Jo 

dTi f dT2 [ dr3H/3(ri)H/3(T2)'H/3(r3) H- 

Jo Jo Jo 

(A31) 


and, using (A31), then 


U{t) = g-itTifi ^[-*/o 


(A32) 


Due that the norm oiUp is very smaller in comparison 
with the norm of 'Hq, a final simplification can be done 
in the case that concerns us. In such case, the value of 


Eq. (A30), for n = 0, applied in the calculation of the 


matrix elements (A25), brings us the result 

ft 


Jo 
= 2 


dTHpir) 




^a.ct' 


Sin 


< 1. 


Observing the integrals in Eq. (A31), we can see that they 


introduce factors that scale the values obtained fro m the 
powers of TList. Because of the small value of Eq. (A33) 


as well as the small values obtained of integrating sev¬ 
eral times expressions like (A33), this scaling produce a 


dynamics very slow. Therefore, if we are not interested 
in the dynamics for very long times, we can neglect in 
Eq. (A31) all the terms with integrals in front of the 


other ones. Finally, the last assumption permit us to 
obtain the leading evolution operator as 


U{t) = e 


— it l-Loc ^—itl-i^ 


(A34) 


3. Relevant evolution operator for the dynamics in 
the decoherence time-scale 


In this section we will obtain the leading evolution 
operator in the decoherence time-scale. In order of that. 


we are using the results of Section A 2 to extract the 
relevant Hamiltonian of such dynamics. The process, to 
obtain that relevant Hamiltonian, consists in splitting 
the Hamiltonian in a strong and a non-commutative 
weak part (i.e. non-commutative with respect to the 
strong part), and then extracting the strong part as 
was developed in Section |A 2[ In that way, with this 
procedure we are defining different adiabatic frames of 


Using the Hamiltonian (|A1|), the form of the exact evo- 

(A35) 


lution operator of the dynamics is 


JJ {t) — 


First, we extract from Eq. (A35) the Zeeman Hamilto¬ 


nian. Thus, we have for the definitions in Eq. (A19) that 
Ha = Hz and Hp =Hd+Hi+Hb- The commutative 
part oi Hp with respect to the Zeeman Hamiltonian 
is Hp =H%+H°i + Hb, where H% = H%^ + H%e = 


H 


0(s 

Da 


(g) -b H°j^"J (g) and H°j = H% + H°j^.^ On 
the other hand, the non-commutative part will heHp = 

{^Da ® + '^De ® ■ 

It is worth that the bath Hamiltonian Hb has eigenval¬ 
ues with absolute values higher that the Zeeman’s ones, 
but that does not affect the approximation due to Hz 
and Hb commute completely between them. 


Now, it is assumed the relationship in Eq. (A33l, where 


{| a)} is in this case the eigenbasis of Hz with eigenvalue 
a. Therefore, in this first adiabatic approximation, the 


evolution operator of the dynamic (A35) can be very 
well written as 

U(t) = (^Da+^De+W/a+^/e+'^^s) (A36) 

(A33) which is well known as the ‘high-held approximation’. 


A second adiabatic frame and approximation can be 
done extra cting the secular dipolar intra-pair Hamilto¬ 
nian from (A36). In this way, we now dehne Ha = TJ^a 
and Hp = + H°^ + H% + Hb- The commuta¬ 

tive part oiHp with respect to H%^ will be written as 
Hp = Hoe + + JJb, where H%^ and H°^ are 

the corresponding commutative part of the dipolar and 
interaction inter-pair Hamiltonians, i.e. and "H/g, 

with respect to the dipolar intra-pair one respec¬ 

tively. It is worth to note that the interaction intra-pair 
and the bath Hamiltonians, i.e. H^^ and Hb, commute 
completely with respect to H'j^a- Accordingly, for the 
non-commutative part we dehne Hp = H%^ + H'j^. 

In order to obtain the form of 'H^g, here we open a paren¬ 
thesis and develop this task in the following. We can see 
from Eqs. (A3), (A4), (A12) and (A13l, that we have to 


consider the problem of truncating the operator 

E E ^0{rA,u-,B,v) 

A;-ii=l,2 B^A;«=1,2 


X T 


(s) 

{2,0}(A,u;B,v)^ 


(A37) 


with regard to 


rp(s) — \ ' rp(s) 
^Da — / , ^(2. 


01 (A)- 


(A38) 


To obtain the effects over the dynamics of truncating, we 
can consider the case where the distance between pairs 
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is enough far to write 

^o{rA,u-,B,v) - ^o{rA,u'-,B,v') = ^oifA.B), ^ U,u',V,v'. 

(A39) 

That approximation as far-away pairs is not the exact 
problem, but to consider the pairs with the actual dis¬ 
tances is a correction to such approximation due that 
only the very-closed pairs do not accomplish exactly the 
relationship (A39|. Otherwise, this approach permits us 


to handle analytically the truncating and to see directly 
the effects over the dynamics, which are very similar to 
the exact calculation. Therefore we can write 


A,B^A 


(A40) 


with 


t(s) 

■{2,0}(A 


- E E T'S 

u=l,2 v=l,2 


0}(A,w;S,i;) ’ 


From the last, we can see that the truncating problem 
is the same that the presented by Keller pS]. where the 
truncated part is 


t(s) 

■(A,B) 


1 

Vi 


2T{i,o}(a)T{i,o}(b) 


+ 2 ('^{i.i}(A)T{i- i}(B) + T{i 


+ 4(T{2,1}(A)T{2 -1}(B) + T{2 -1}(A)T{2,1}(B)) 




where 


-i(s) 

-{1,0}(A) 


t(®) 


r(^) I t(^) 

v.(A,l) z(A.2)> 


(A41) 

(A42a) 


OjO 0/0 
^De-> ^Dq 


= 0 . 


with the commutation property 
It is worth to mention that, in the analysis presented in 
this work, we are not concerned for the explicit form of 

ttl^. 

As in the first case, we will assume the relationship in 
Eq. (A33), where {| a)} is in this case the eigenbasis of 
with eigenvalue a. 

Finalizing this second adiabatic approximation, we con¬ 
clude that the evolution operator of the dynamics can be 
written as 

U{t) ^ e-*‘(«De+«?a+w?.+WB)_ (^46) 

Following the extraction of adiabatic frames from the 
dynamics, we can perform a third approximation, and a 
new frame definition, extracting from the evolution op¬ 
erator (A46) the truncated dipolar inter-pair Hamilto¬ 
nian In order of that, we define Ha = H%e 

Hp = H'j^ -I- H/g + Hb- For this case, the commuta¬ 
tive part oi Hp with respect to H'j^^ will be written as 
Up = TL'jg^+'Hll+HB, where and are the com¬ 
mutative parts of the interaction intra-pair and the once- 
truncated interaction inter-pair Hamiltonians, i.e. H^^ 
and H/g, with respect to the truncated dipolar inter¬ 
pair Hamiltonian respectively. Besides, we have 

that H^g commute with Hb- Accordingly, for the non- 
commutative part we define Hp = H'j^^ -t- "Hjg. 

Here we open a new parenthesis to obtai n, in t he fo llow- 
ing, the form of H^^. From Eqs. (A6al and (A14), we 


can write the intra-pair interaction Hamiltonian as 


'TJV — 

^la = ~ 


3V6 0o(d) 


Et 


is) 

k 


-Ab) 


(A47) 


where are defined 


{1,±1}(A) - + ^^±iA) - + ^ V ±(A1) + ^±iA 


is) _ 


1 / 


is) 


. 2 )) ’ 

(A42b) 


and the operators are defined in Eq. ( All|. 

To close this parenthesis, using the ideas exposed in the 
last paragraphs, we now are able to write the expression 
for 77^)6 ■ Using the following definitions 


= Ho (a) ET 


(«) 

De’ 


(A43) 


with Et^^ ^ Ea T 


De(A) 


, and 


As) _ 

-De{A) 


-E 

B^A 


a 


rA,B 


As) 

■ (A,B)^ 


where a is the shortest distance between pairs, we obtain 
therefore 


H° 


V6. 


De — 2 Uo(a) ET^g 0 1 


ib) 


(A45) 


-\(s) 


-Et 


(s) 

{2.01(A) ' 


^ — ikAa 




E[^^Vbg+.MbI? 


(A48a) 


(A48b) 


Therefore, the last truncating problem is equival ent to 
obtain the truncated part of each Tg in Eq. (A47) 
with respect to H^g. If it is defined a maximum value 
of k, namely k, such that |fcAa| <C 1 and therefore 

^—ikAa 


(A44) 0. 


■^{s) ^ 
-k<k 


= E 


1, V fc < fc, we can verify from Eq. (|A48a|) that 

is) 


T' 

A {2.01(A) 


= and thus 


4/0 fy(^) 

^k<k 


Finally, to close the last pare nthes is, using the last con¬ 
clusions we can neglect in Eq. (A47) the terms with k > k 
and then define 


-l/U — 

^la = ~ 


3V6Ho(d) 


Et 


k 


-Ab) 


(A49) 


k<k 
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It is worth to mention that, as we will see in the follow¬ 
ing, we are not concerned for the explicit form of ’Hj®. 
Again, we will neglect the evolution o perato r produced by 
Up assuming the relationship in Eq. (A33), where {la)} 
is in this case the eigenbasis of eigenvalue a. 

Finalizing this third adiabatic approximation, we con¬ 
clude that the evolution operator of the dynamics can be 
written as 


U{t) ^ e +«i3) ^ 


Now we can assume the following relationship between 
thus the decoherence produced 


norms, 


'0 




< 


njO 

^la 


by the Hamiltonian will be faster than the produced 


by 'Hjg. Accordingly, for a shake of simplicity in the 
analytical calculation of the dynamics, we will neglect 
the dynamics generated by 

Therefore, to conclude this section, we define the rele¬ 
vant evolution operator, which governs the dynamics in 
the decoherence time-scale, as the adiabatic expression 




(A51) 


where Hz, H^,^ and H 
with 7 }°^ and Hb, but 


commute between them and 




^0. 
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